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We consider the nature of the fluid-solid phase transition in a polydisperse mixture of hard 
spheres. For a sufficiently polydisperse mixture (a > 0.085) crystallisation occurs with simultaneous 
fractionation. At the fluid-solid boundary, a broad fluid diameter distribution is split into a number 
of narrower fractions, each of which then crystallise. The number of crystalline phases increases with 
the overall level of polydispersity. At high densities, freezing is followed by a sequence of demixing 
transitions in the polydisperse crystal. 
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' Equal-sized hard spheres constitute probably the simplest example of a purely entropic material. In a hard-sphere 
system there is no contribution to the internal energy U from interparticle forces so that U is a constant, at a fixed 
temperature. Minimising the free energy, F = U — TS, is thus simply equivalent to maximising the entropy S. Conse- 
' quently, the structure and phase behaviour of hard spheres is determined solely by eaitropy. Although the hard-sphere 
. model was originally introduced as a mathematically simple model of atomic liquidsEl recent work has demonstrated it 
7— I ' usefulness as a basic model for complex fluids! Suspensions of submicron poly(methyl methacrylate) or silica colloids, 
coated with a thin polymeric layer so that strong repulsions dominate the attractive dispersion forces between the 
colloidal cores, behave in many ways as hard spheres. In the last decade, a great deal of effort has been devoted 
to systematic studies of such colloidal 'hard-sphere' systems! For sufficiently monodispersej-colloids, crystallisation 
is observed at densities similar to those predicted by computer simulation for hard spherearEl. Measurements gf the 
osmotic pressure and compressibility show a similar dramatic agreement with predicted hard sphere propertiesH. 

There is however one important and unavoidable difference between colloids and the classical hard sphere model 
which is frequently overlooked. Whereas the spheres in the classical model are identically-sized (i.e. monodisperse) 
colloidal particles have an inevitable spread of sizes which is most conveniently characterised by the polydispersity, 
I 1 <7, defined as 
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where 

pW = J dRp(R)R n (2) 

with p(R) the density distribution and p = J dRp(R). 

Recent work has revealed that as soon as a hard sphere suspension is allowed to enjoy a significant degree of poly- 
dispersity, several interesting new phenomena arise. Experiments find that the crystallisation transition is suppressed 
altogether at a w 0.12,while samples with a = 0.075 crystallise only slowly in the coexistence region and not above 
the melting transitional. Other exampl£s_of the effects of polydispersity are the appearance of a demixing transition 
in a polydisjperse fluid of hard spheresau and the observation of a liquid-vapour transition in polydisperse adhesive 
hard spheresa. 

The effect— pf polydispersity on— the crystallisation of hard . sph ere colloids has been investigated by computer 
simulationcTO, density functionallla't^l, and analytical theoriealjlla. The picture that emerges is remarkably con- 
sistent. All calculations find that the fluid-solid phase transition vanishes for polydispecsities above a certain critical 
level a c . The phase diagram for small polydispersities (a < a c ) has been rationalisedtJ in terms of the appearance 
of an additional high density crystal-to-fluid transition, in a polydisperse system. While, at low polydispersities 
hard spheres display, with increasing density, the conventional fluid-to-crystal transition, at higher polydispersities 
re-entrant behaviour is predicted. The two freezing transitions converge to a single point-iri the (</>, a) plane which 
is a polydisperse analogue of the point of equal concentration found in molecular mixturesEJj. At this singularity, the 
free energies of the polydisperse fluid and crystal phases are equal. 
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The purpose of this note is to examine the fate of a highly polydisperse (er > a c ) hard sphere fluid. Previous 
theoretical research has not been able to identify a fluid-solid transition for a > a c so it is generally believed that 
the equilibrium phase is disordered at all densities up to close packing. Several years ago, Puseyt-3 suggested that a 
highly polydisperse suspension might crystallise by splitting the broad overall distribution into a number of narrower 
distributions of polydispersity <7 S , each of which could be accommodated within a single crystalline phase. Each 
crystal would then have a correspondingly different mean size with the number of crystalline phases increasing with 
the overall polydispersity. For fractionated crystallisation to occur the total free energy of the set of multiple crystals 
must be lower than that of the equivalent polydisperse fluid. This can only happen if the reduction in free energy 
as particles are removed from a fluid and placed in a crystal is sufficiently large to exceed the loss of entropy of 
mixing as the distribution is partitioned. This is a delicate balance and it is far from obvious where the result lies. 
Pusey, for instance, originally suggested that fractionation would generate crystals with a polydispersity of a c so as 
to minimise the number of crystal phases required and the subsequent loss of entropy of mixing. However, as noted 
above, at a c the free energy of the polydisperse fluid and crystal phases are equaltJ and so there is no driving force 
for fractionation.!—. 

In earlier worked the possibility of fractionated crystallisation for a > a c was considered but no conditions where 
two crystal phases could coexist could be found. Here we re-examine the stability of a polydisperse hard-sphere fluid 
using a much simpler approach. Rather than solving the equations of phase equilibria in a polydisperse system we 
restrict ourselves to the easier task of comparing the free energies, at the same density and temperature, of crystal 
and fluid phases. We find, in agreement with PuseyEZI, that fractionation occurs in polydisperse hard sphere mixtures 
but that the polydispersity of the resulting crystals is substantially less than a c . 

The rest of the paper is organised as follows: in section || we present our model for the free energies of the 
polydisperse fluid and crystal phases. The stability diagram is present in section III. Finally, in section IV we 
summarise and conclude. 



II. THE MODEL 



Our model consists of N hard sphere particles in a volume V, at an overall density of p — N/V. Each particle has 
a diameter R drawn from some overall distribution p(R) so that p = J dRp(R). Previous work has suggestedE-3 that 
the thermodynamic properties of a polydisperse system are relatively insensitive to the detailed shape assumed for 
the diameter distribution, at least when the polydispersity is small (a < 0.1). Indeed the phase behaviour has been 
foundEj, to a rather good approximation, to be a function only of the first three moments of the diameter distribution, 
p, R, and the normalised width a. For a broad diameter distribution, higher moments will presumably also need to 
be considered. Consequently, the phase behaviour will become increasingly sensitive to the shape of the distribution 
assumed. However we expect that the width of the distribution will still play an important if not the dominant role 
in determining the generic features of the phase behaviour of a polydisperse system. Since our concern here is with 
establishing these general features of the polydisperse phase diagram we have taken p(R) as a simple rectangular 
distribution, 

( R<R(1- w/2) 

p(R)=pt (Rw)- 1 R(l -to/2) < R < R(l + w/2) (3) 
[O R>R(l+w/2) 

where R is the mean diameter and the normalised width w is related directly to the conventional polydispersity a by 
w = 2\/?><7. It will prove convenient to choose as variables the mean diameter R, the polydispersity a and the volume 
fraction = (vr/6) / dRp(R)R? = {ir/&)pR?{l + w 2 /4). 

Our first task is to calculate the thermodynamics of the fluid state. The free energy density (/ = F/V) may be 
conveniently split into ideal and excess portions 

./' f d I ./" ' . (4) 

The ideal part of the free energy density is given by a straightforward generalisation of the free energy of an ideal gas 
as 

0f id = J dRp(R)(lnp(R)-l) (5) 
with j3 — 1/(&bT). For the case of the rectangular distribution, the ideal free energy density is simply 
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f3r d = p(lnp-l)-pHRw). (6) 

We calculate the excess free energy density, f ex , of the polydisperse fluid from the general equation of state given 
by Salacuse and StellEl This equation of state is a straightforward generalisation of the highly successful EOS for 
an arbitrary hard sphere mixtureciJ first considered by Mansoori et al. The Mansoori EOS has been checked against 
simulation data by a number of authors and agreement is generally excellent. The free energy per particle (J- — F/N) 
is then Tj = {ff + If )/p. 

Now we consider the free energy of the fractionated solid phase. Since as mentioned above, we expect the number 
of crystals to depend upon the overall polydispersity a of the parent distribution, we consider the general case of P 
coexisting crystals. In order to establish some relation between the parent distribution p(R) and the sub-distribution, 
Pi{R), found in the i th crystal we make two assumptions. First we choose the distribution in each crystal to be rectan- 
gular and second, to keep the theory computationally manageable, we force the width (or polydispersity) to be equal in 
each of the P coexisting phases. With these two restrictions the housekeeping of the fractionation process is straight- 
forward. The i th crystal, for instance, contains all iVj particles whose diameters lie between R [l — w s (l + -j — i)] 
and R [l — w s (4p — i)] with w s , the width of the crystal distribution, given by w s = w/P. The mean diameter of 
particles in the i th crystal is then simply given as 



Ri — R 
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l__i(l + P _2i) (7) 
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and the polydispersity of each crystal is a s = a/ P. The diameter distribution in the i th crystal is then 

R<Ri(l-w„/2) 

p i {R) = p i { (RtWs)- 1 Ri(l-w s /2)<R<Ri(l + w s /2) (8) 
R > R t (l + w s /2) 

at an overall density of pi = Ni/Vi with Vi the volume of the i th phase. Although the number of particles in each 
crystal phase is the same for all crystals and equal to N/P (since the polydispersity of each crystal is assumed equal) 
the volume Vi of each phase is so far undetermined. To fix the volume, or equivalently the number density p^, of 
each crystal we use the fact that the crystals are in equilibrium with each other. This constraint is straightforward 
to apply since, in our model of the partition, no particle may exist in anymore than one crystal. Consequently we do 
not need to consider exchange of particles. Equilibrium simply requires that the osmotic pressure of each crystal, IT, 
be fixed and equal to an external pressure of, say, II 

TU(jH,ir. i R i )=Il. (9) 
With the density of each crystal phase determined, the total number density of the set of P coexisting crystals is then 

Ps = — p— — (10) 

and the corresponding composite packing fraction is = (7r/6)p s _R 3 (l + u> 2 /4). The mean free energy per particle is 

where fi is the free energy density of the i th crystal. 

To complete our model for the fractionated phases we need an expression for the excess free energy of a polydisperse 
crystal. For that purpose we use a scaled particle theory. Since the details have been described elsewhereEJ we just 
outline the approach here. The model exploits the picture of particles in a solid as being confined in cells or cages 
formed by the neighbours from which they can not escape. While, in a monodisperse crystal the cells are all identical, 
in a polydisperse crystal the size and shape of the cell varies. The idea is to recognise that one can allow for the 
variability in the cell shape, to a first approximation, by considering just a finite number of different-sized neighbours. 
The number of different-sized particles is determined by the accuracy with which the properties of the finite mixture 
approximate those of the polydisperse system. If the excess free energy of the polydisperse system is a function 
of the first p moments of the diameter distribution [p/2] different-sized spheres are needed to match the diameter 
moments of the continuous distribution ([x] is the smallest integer not less than x). For a polydisperse system of 
hard spheres, highly successful rpean- field theories such as the Percus-Yevick approximation suggest that only four 
diameter moments are significant. In this case, the properties of a polydisperse system should be well approximated 
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by a binary mixture of spheres, chosen so that the first four moments of the polydisperse and binary distributions 
are equal. The two diameter distributions are 'equivalent' and all excess properties are, to a first approximation, 
equal. For the rectangular distribution the equivalent binary distribution has an equal number of large and small 
spheres with diameters Ri(l + cr s ) and Ri(l — er s ) respectively. Specifically then, wc equate the excess free energy of 
the polydisperse crystal with that of an equivalent binary substitutional crystal__For the evaluation of the latter we 
take advantage of the analytical expressions quoted by Kranendonk and FrenkelEil as fits to Monte Carlo simulation 
results. Combining the excess free energy with the ideal free energy density, 

0ft d = Pi {\n Pi -I) - Pi HR lWs ) (12) 

yields the free energy density fa of the i th crystal. Finally, we note that our model for the polydisperse crystal 
implicitly assumes that all different-sized spheres are placed randomly on the sites of a common fee lattice. At high 
polydispersities, the smallest particles in the distribution could also be accommodated within the interstices of a 
crystal of larger particles. The current calculations ignore this possibility and so are valid only for polydispersities of 
a < 0.3. 

Putting everything together, we now calculate the difference free energy per particle, AT P = T p — Ff, between 
the polydisperse fluid and the fractionated set of P crystals as follows. The overall particle density distribution of 
either the set of crystals or the polydisperse fluid is characterised by identical values for the overall volume fraction 
(j>, the polydispersity a and the mean diameter R. For fixed <f> = <f>, a and R, the free energy of the fluid Tj is 
calculated from Eqs. ^ and the Mansoori EOS. Combining Eqs. || and |l^ we see that the overall volume fraction <f> 
of the solid phases (or equivalently the mean density p s ) may be considered as a function of the osmotic pressure 
IT. Inverting this relation yields the osmotic pressure II and so, from Eq. [)| the densities pi of each of the coexisting 
crystals which taken together have an overall volume fraction of <fi- Knowing the density, polydispersity a s = cr/P 
and mean diameter Ri (from Eq. [7]) of each crystal it is straightforward to calculate the corresponding free energy 
density fi and the mean free energy T p of the multiple crystals from Eq. nil 



III. RESULTS 



We now present in detail results for the stability of a system of polydisperse hard spheres obtained from the theory 
described above. In Fig. [j], we show the free energy difference, /3AJF P , between the fractionated crystals and the fluid 
phase as a function of polydispersity, for the representative density of <p = 0.60. The first question to be addressed is 
the relative stability of the various fractionated crystals and the fluid phase. 

Referring to Fig. EL we see that for low levels of polydispersity, at this density, the crystal has a significantly lower free 
energy than the fluid phase. However, with increasing polydispersity, the stability of the polydisperse crystal reduces 
rapidly until for a > 0.084 the unfractionated polydisperse crystal is unstable relative to the fluid phase. The origin 
of this behaviour lies in the different ways polydispersity affects the maximum packing fraction, or equivalently the 
free volume, of ordered and disordered structures. In a fluid, smaller particles are free to pack in the cavities between 
larger particles and so the free volume and thus the entropy increases with polydispersity. Conversely, the periodicity 
of a crystal causes the free volume and entropy of a polydisperse crystal to decrease with increasing polydispersity. 
Consequently the free energy of a fixed-density polydisperse crystal will diverge at a critical level of polydispersity at 
which the fluid free energy will remain finite. Second, we note from Fig. EJ that at low polydispersities, fractionating 
the distribution into two or more crystals always raises the free energy of the solid state. This is because of the loss 
of entropy of mixing as the diameter distribution is split up. The contribution to the free energy difference, (3AT P , 
due to the loss of entropy of mixing is simply: 

(3AT p lx = \nP (13) 

for fractionation into P crystal phases which is very close to the differences seen in Fig. [I] as a — > 0. However at finite 
polydispersities, Eq. does not give a reasonable estimate of the free energy differences between the various sets of 
crystals. As is evident, the free energy j3AT p reduces with increasing polydispersity. The reduction being smaller 
as the number of crystals considered increases. This is because the polydispersity of each fractionated crystal drops 
as the total number of crystals, between which the overall distribution is split, increases. For instance, fractionating 
a distribution with a polydispersity of 6% between two phases would result in crystals with 3% polydispersity, three 
phases would lead to 2% polydispersity and so on. As the polydispersity of each of the individual crystals decreases 
the influence of the divergence of the free energy at the close-packing limit is reduced and so the free energy is less 
affected by the overall polydispersity. As a consequence the fractionated two-crystal phase takes the place of the 
polydisperse single-phase crystal, as the state of lowest free energy, for a > 0.078 and remains the most stable phase 
until a = 0.099, at which point the fluid appears. 
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By a procedure similar to that described above, we have mapped out the range of stability of the polydisperse hard 
sphere system for a < 0.20. Our results are shown in Fig. |2[ The regions are labelled by the number of coexisting 
crystals which minimise the total free energy of the polydisperse system while the boundaries indicate the specific 
densities and polydispersities at which competing phases have equal free energies. As is evident, there are appreciable 
areas of density and polydispersity where fractionation into multiple crystals is predicted. 

At low volume fractions, <f> < 0.516, the polydisperse fluid is the most stable phase at all polydispersities. However 
with an increase in packing fraction a crystalline phase takes the place of the polydisperse fluid as the most stable 
phase. The nature of the crystalline phase depends on the magnitude of the polydispersity a. So while for a < 0.085 
we have crystallisation into the usual single polydisperse crystal, with increasing polydispersity the fluid fractionates 
upon solidification into a rapidly increasing number of crystals. The diagram gets more and more complicated as 
the polydispersity grows since, by necessity, the number of crystalline phases required to accommodate the increasing 
width of the diameter distribution increase. So while for polydispersities in the range 0.085 < a < 0.117 two 
fractionated crystals are first formed, three crystals are required for 0.117 < a < 0.145, four for 0.145 < a < 0.171, 
five for 0.171 < a < 0.195 while at a — 0.20 complete crystallisation requires six different crystal phases. Furthermore, 
we note that as the polydispersity increases the packing fraction at which the fluid becomes unstable increases from 
<\> = 0.516 at a — to <p — 0.665 at a = 0.20. The degree of fractionation predicted varies with the overall level of 
polydispersity. Referring to Fig. ||, we see that the polydispersity of the crystal(s) at the fluid-solid boundary reaches 
a maximum of 0.085 before on average falling with increasing fluid polydispersity, although in a far from continuous 
fashion. Where fractionation occurs, Fig. ||| reveals that the polydispersity of each crystal formed is considerably lower 
than the critical polydispersity, estimated here as a c ~ 0.085. 

In addition to the fluid- multiple crystal transitions discussed above, Fig. || reveals that there is a further sequence 
of demixing transitions in the solid phase. Compressing a polydisperse crystal causes a solid-state phase separation 
in which the diameter distribution is fractionated. For instance a polydisperse crystal with a = 0.08 is stable up to a 
density of </> = 0.597. At this point the system separates into two crystals, each of polydispersity a = 0.04, which are 
stable until a total density of <f> — 0.667 is reached. Further compression results in three crystals, followed at still high 
densities by a cascade of demixing transitions. The explanation for this behaviour lies in the effect of polydispersity 
on the limiting packing fraction of a crystal. As mentioned previously, the density at which the free energy of a 
polydisperse crystal starts to diverge reduces with increasing polydispersity. As a polydisperse crystal is compressed 
there comes a point at which the reduction in the excess free energy which occurs when the polydispersity is reduced 
is sufficient to exceed the increase in the ideal free energy of mixing after fractionation. At this point the crystal 
phase separates. 

Finally, we emphasise that we have calculated only the stability boundaries, not the full equilibrium phase diagram 
of a system of polydisperse spheres. The reason for this is that a full calculation of the phase equilibria, particularly 
in view of the large number of competing phases found in this work, would be a large and complicated problem. 
However we expect our stability diagram (Fig. ^|) to be a useful guide to the form of the phase diagram. Confirmation 
of this point of view is provided by the comparison.-seen in Fig. |J, between the current model and the fluid-single 
crystal coexistence region established in earlier workli3. The dashed curves indicate the positions of the polydisperse 
cloud-point boundaries and were calculated by a novel moment projection method using, as input, polydisperse free 
energies similar to those described above. Referring to Fig. fl|, we see that the stability boundary lies approximately 
midway between the two coexistence densities. Inspection of Fig. reveals, rather intriguingly, that the polydisperse 
point of equal concentration, marked by the large circle in Fig. |J is practically coincident with the lowest density 
at which the binary system of fractionated crystals is stable. This seems to be a chance occurrence.. However, the 
proximity of the two points suggests that, at equilibrium, the re-entrant melting predicted in Ref.E3 will in fact be 
pre-empted by a solid state phase separation. 

IV. DISCUSSION AND CONCLUSIONS 

We have analysed the stability of a polydisperse fluid of hard spheres with respect to a process of simultaneous 
fractionation and crystallisation. Although the model is quite simple we find several interesting features. In particular, 
we predict that for a sufficiently polydisperse system of hard spheres fractionation occurs upon solidification. The 
fluid diameter distribution separates into a number of fractions of narrower polydispersity which then crystallise. 
Furthermore, we find that compressing a polydisperse crystal induces a sequence of demixing transitions in the 
crystal. 

To keep the theory presented here manageable we have been forced to make a number of assumptions. We have, 
for instance, imposed a somewhat arbitrary model for the fractionation process. In particular we have assumed that 
a rectangular diameter distribution fractionates into a number of equal-width daughter distributions, each of which 
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is also rectangular in form. In principle, we could relax the constraint of equal width. However, since the densities 
of the coexisting crystals are similar it seems likely that this modification would have little effect. The more critical 
restriction is probably the assumption made for the form of the daughter distribution. However, in view of the relative 
insensitivity of polydisperse systems to the detailed form chosen for the diameter distribution, we expect that the 
current predictions will be in at least qualitative if not quantitative agreement with more rigorous calculations, when 
available. 

We now speculate briefly on the feasibility of observing solid-state fractionation in experiments on hard-sphere 
colloids. The first point to note is that fractionation requires colloid diffusion over distances comparable to the size 
of the growing crystallite. The rate of such large scale diffusive motion reduces with increasing colloid density, as a 
result of caging effects, and essentially vanishes at the glass transition density, 4> g a. For uniformly-sized hard spheres 
simulationsEa find a long-lived metastable glassy state around (j) g = 0.58 which persists up to the random close packing 
limit at 4> rc p = 0.64. The glass transition in a-system of polydisperse hard spheres has not been studied to the best of 
our knowledge although Schaertl and Silcscua have estimated <fi rcp as a function of polydispersity from simulation. 
Assuming a simple linear relationship between the two densities gives a crude estimate of the glass transition density 
as 4> g {<j) ~ (0.58/0.64)</> rC p(er). Referring to Fig. |], we see that if this estimate for the glass transition is used, then 
fractionated crystallisation might be experimentally observable for polydispersities in the narrow range 0.08 < a < 0.11 
where diffusive motion is still possible although slow. If conversely, the dynamics of the fractionation process turn out 
to be appreciably slower than the experimental timeframe then experiments should follow our earlier predictionst3 
and exhibit a polydisperse point of equal concentration. A clear resolution of this ambiguity must await experiment. 

Finally, our results suggest a possible explanation for the apparently anomalous simulations reported by Bolhuis and 
Kofkd_J. In contrast with other theoretical work, this study found that: (i) the particle size distributions in the fluid 
and crystal phases were significantly different and (ii) that the coexistence region although initially narrowing with 
increasing polydispersity finally widened. The phase diagram was traced out by integrating along the fluid-crystal 
coexistence line. Implicit in this approach is the assumption that only one crystalline phase exists at the fluid-solid 
boundary. Referring to Fig. @, we envisage that at the point where two polydisperse crystals become stable the 
simulation will follow one of the two fluid-solid branches until it halts at the discontinuity marking the start of the 
three-crystal region. With this picture in mind, we may readily account for the observation that at the point, where 
the simulation could no longer follow the fluid-solid boundary, the crystal polydispersity was almost exactly half that 
of the fluid value. At the same time, we note that their limiting polydispersity (cr = 0.118) lies close to our estimate 
for the highest polydispersity (0.117) at which two polydisperse crystals remain the equilibrium phase. 
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Bartlett, Fig. 1 




FIG. 1. The free energy difference per particle, at a constant packing fraction of <f> = 0.6, between polydisperse crystalline 
and fluid phases as a function of polydispersity. The number of coexisting crystals is detailed in the legend. The arrows mark 
the polydispersity limits between which the fractionated binary crystals are stable. 
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Bartlett, Fig. 2 




FIG. 2. The stability diagram of a polydisperse hard sphere mixture. In each region of volume fraction <j> and polydispersity 
a the phase with the lowest free energy is indicated. Crystalline regions are labelled by the number of coexisting crystal phases. 
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Bartlett, Fig. 3 




FIG. 3. The polydispersity, a s , of the fractionated crystalline phase formed at the fluid-solid boundary as a function of the 
fluid polydispersity a. 
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Bartlett, Fig. 4 




FIG. 4. A comparison between the current stability diagram (solid lines) and the fluid-solid coexistence predicted in 
(dashed lines). The large circle marks the position of the polydisperse point of equal concentration. The thick dotted line 
denotes the estimated location of the polydisperse glass transition. 
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